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Spin pumping is the emission of a spin current by a magnetization dynamics while spin transfer 
stands for the excitation of magnetization by spin currents. Using Onsager's reciprocity relations 
we prove that spin pumping and spin-transfer torques are two fundamentally equivalent dynamic 
processes in magnetic structures with itinerant electrons. We review the theory of the coupled 
motion of the magnetization order parameter and electron for textured bulk ferromagnets (e.g. 
containing domain walls) and heterostructures (such as spin valves). We present first-principles 
calculations for the material-dependent damping parameters of magnetic alloys. Theoretical and 
experimental results agree in general well. 
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I. INTRODUCTION 
A. Technology Pull and Physics Push 

The interaction between electric currents and the magnetic order parameter in conducting magnetic micro- and 
nanostructures has developed into a major subfield in magnetism^. The main reason is the technological potential 
of magnetic devices based on transition metals and their alloys that ope rate at ambient temperatures. Examples 
are current-induced tunable microwave generators (spin-torque oscillators) 23 , and non-volatile magnetic electronic 
architectures that can be randomly read, written or programmed by current pulses in a scalable manner 4 . The 
interaction between currents and magnetization can also cause undesirable effects such as enhanced magnetic noise 
in read heads made from magnetic multilayers-. While most research has been carried out on metallic structures, 
current-induced magnetization dynamics in semiconductors^ or even insulators^ has been pursued as well. 

Physicists have been attracted in large numbers to these issues because on top of the practical aspects the underlying 
phenomena are so fascinating. Berger^ and SlonczewskP are in general acknowledged to have started the whole field 
by introducing the concept of current-induced magnetization dynamics by the transfer of spin. The importance o f their 
work was fully appreciated only after experimental confirmation of the predictions in multi-layered structure d 10 * 1 H 
The reciprocal effect, i .e.. the generation of currents by magnetization dyn amics now called spin pumping, has 
been expected long age) 12 * 13 1 , but it took some time before Tserkovnyak et a/! 14 * 15 1 developed a rigorous theory of 
spin-pumping for magnetic multi-layers, including the associated increased magnetization dampin g 16 " 18 !. 

B. Discrete versus Homogeneous 

Spin-transfer torque and spin pumping in magnetic metalli c multi- layers are by now relatively well understood 
and the topic has been covered by a number of review articles^ 1121221 . It can be understood very well in terms of 
a time-dependent extension of magneto-electronic circuit theory^ 2 ^, which corresponds to the assumption of spin 
diffusion in the bulk and quantum me chanical boundary conditions at interfaces. Random matrix theory^ can be 
shown to be equivalent to circuit theor y 19 * 23 * 2 [ The technologically important current-induced switching in magnetic 
tunnel junctions has recently been the focus of attention*^. Tunnel junctions limit the transport such that circuit 
issues are less important, whereas the quantum-mechanical nature of the tunneling process becomes essential. We 
will not review this issue in more detail here. 

The interaction of currents and magnetization in continuous magnetization textures has also attracted much interest, 
partly due to possible applications such as nonvolatile shift registers^. From a formal point of view the physics 
of current-magnetization interaction in a continuum poses new challenges as compared to heterostructures with 
atomically sharp interfaces. In magnetic textures such as magnetic domain walls, currents interact over length scales 
corresponding to the wall widths that are usually much longer than even the transport mean-free path. Issues of the 
in-plane vs. magnetic- field like torque^ and the spin- motive force in moving magnetization textures^ took some time 
to get sorted out, but the understanding of the complications associated with continuous textures has matured by now. 
There is now general cons ensus about the physics of current-induced magnetization excitations and magnetization 
dynamics induced current ^ 29 * 30 *. Nevertheless, the similarities and differences of spin torque and spin pumping in 
discrete and continuous magnetic systems has to our knowledge never been discussed in a coherent fashion. It has 
also only recently been realized that both phenomena are directly related, since they reflect identical microscopic 
correlations according to the Onsager reciprocity relationship. 

C. This Chapter 

In this Chapter, we (i) review the basic understandings of spin transfer torque vs. spin pumping and (ii) knit together 
our understanding of both concepts for heterogeneous and homogeneous systems. We discuss the general phenomenol- 
ogy guided by Onsager's reciprocity in the linear response regime^. We will compare the in- and out-of-plane spin 
transfer torques at interfaces as governed by the real and imaginary part of the so-called spin-mixing conductances 
with that in textures, which are usually associated with the adiabatic torque and its dissipative correction^, usually 
described by a dimensionless factor /3 in order to stress the relation with the Gilbert damping constant a. We argue 
that the spin pumping phenomenon at interfaces between magnets and conductors is identical to the spin-motive force 
due to magnetization texture dynamics such as moving domain walls^. We emphasize that spin pumping is on a 
microscopic level identical to the spin transfer torque, thus arriving at a significantly simplified conceptual picture of 
the coupling between currents and magnetization. We also point out that we are not limited to a phenomenological 



3 



description relying on fitting parameters by demonstrating that the material dependence of crucial parameters such 
as a and /3 can be computed from first principles. 

II. PHENOMENOLOGY 

In this Section we explain the basics physics of spin-pumping and spin-transfer torques, introduce the dependence 
on material and externally applied parameters, and prove their equivalence in terms of Onsager's reciprocity theorem. 

A. Mechanics 

On a microscopic level electrons behave as wave-like Fermions with quantized intrinsic angular momentum. However, 
in order to understand the electron wave packets at the Fermi energy in high-density metals and the collective motion 
of a large number of spins at not too low temperatures classical analogues can be useful. 

Spin transfer torque and spin pumping are on a fundamental level mechanical phenomena that can be compared 
with the game of billiards, which is all about the transfer of linear and angular momenta between the balls and 
cushions. A skilled player can use the cue to transfer velocity and spin to the billiard ball in a controlled way. The 
path of the spinning ball is governed by the interaction with the reservoirs of linear and angular momentum (the 
cushions and the felt/baize) and with other balls during collisions. A ball that for instance hits the cushion at normal 
angle with top or bottom spin will reverse its rotation and translation, thereby transferring twice its linear and angular 
moment to the frame of the billiard. 

Since the work by BarnettP^and Einstein-de Haas^ almost a century ago, we know that magnetism is caused by the 
magnetic moment of the electron, which is intimately related with its mechanical angular momentum. How angular 
momentum transfer occurs between electrons in magnetic structures can be imagined mechanically: just replace the 
billiard balls by spin polarized electrons and the cushion by a ferromagnet. Good metallic interfaces correspond to a 
cushion with high friction. The billiard ball reverses angular and linear momentum, whereas the electron is reflected 
with a spin flip. While the cushion and the billiard table absorb the angular momentum, the magnetization absorbs 
the spin angular momentum. The absorbed spins correspond to a torque that, if exceeding a critical value, will set 
the magnetization into motion. Analogously, a time-dependent magnetization injects net angular momentum into 
a normal metal contact. This "spin pumping" effect, i.e. the main topic of this chapter, can be also visualized 
mechanically: a billiard ball without spin will pick up angular momentum under reflection if the cushion is rotating 
along its axis. 

B. Spin-transfer Torque and Spin-pumping 

Ferromagnets do not easily change the modulus of the magnetization vector due to large exchange energy costs. 
The low-energy excitations, so-called spin waves or magnons, only modulate the magnetization direction with respect 
to the equilibrium magnetization configuration. In this regime the magnetization dynamics of ferromagnets can be 
described by the Landau-Lifshitz-Gilbert (LLG) equation, 

rh = —7m x H c ff + am x rh, (1) 

where m(r,i) is a unit vector along the magnetization direction, rh = dm/dt, 7 = g*/!^/^ > is (minus) the 
gyro-magnetic ratio in terms of the effective g-factor and the Bohr magneton and a is the Gilbert damping 
tensor that determines the magnetization dissipation rate. Under isothermal conditions the effective magnetic field 
H e ff = — 5F [m] /5(M s m) is governed by the magnetic free energy F and M s is the saturation magnetization. We will 
consider both spatially homogeneous and inhomogeneous situations. In the former case, the magnetization is constant 
in space (macrospin), while the torques are applied at the interfaces. In the latter case, the effective magnetic field 
H c ff also includes a second order spatial gradient arising from the (exchange) rigidity of the magnetization and torques 
as well as motive forces that are distributed in the ferromagnet. 

Eq. (JlJ can be rewritten in the form of the Landau-Lifshitz (LL) equation: 

(l + a 2 ) m = —7m x H c ff — 7am x (m x H cff ) . (2) 

Additional torques due to the coupling between currents and magnetization dynamics should be added to the right- 
hand side of the LLG or LL equation, but some care should be exercised in order to keep track of dissipation in 
a consistent manner. In our approach the spin-pumping and spin-transfer torque contributions are most naturally 
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added to the LLG equation Q, but we will also make contact with the LL equation ([2]) while exploring the Onsager 
reciprocity relations. 

In the remaining part of this section we describe the extensions of the LLG equation due to spin-transfer and 
spin-pumping torques for discrete and bulk systems in Sec. |II B l| and Sec. |IIB 2[ respectively. In the next section we 
demonstrate in more detail how spin-pumping and spin-transfer torque are related by Onsager reciprocity relations 
for both discrete and continuous systems. 



1. Discrete Systems 



Berger and Slonczewski predicted that in spin-valve structures with current perpendicular to the interface planes 
(CPP) a dc current ca n ex cite and even reverse the reverse the relative magnetization of magnetic layers se parated by 
a normal metal spaceiP^. The existence of this phenomenon has been amply confirmed by experiment d 10 * 11 * 20 * 37 * ^. 
We can understand current-induced magnetization dynamics from first principles in terms of the coupling of spin- 
dependent transport with the magnetization. In a ferromagnetic metal majority and minority electron spins have often 
very different electronic structures. Spins that are polarized non-collinear with respect to the magnetization direction 
are not eigenstates of the ferromagnet, but can be described as a coherent linear combination of majority and minority 
electron spins at the given energy shell. If injected at an interface, these states precess on time and length scales that 
depend on the orbital part of the wave function. In high electron-density transition metal ferromagnets like Co, Ni, 
and Fe a large number of wave vectors are available at the Fermi energy. A transverse spin current injected from a 
diffuse reservoir generates a large number of wave functions oscillating with different wave length that lead to efficient 
destructive interference or decoherence of the spin momentum. Beyond a transverse magnetic coherence length, which 
in these materials is of the order of the Fermi wave length, typically around 1 nm, a transversely polarized spin current 
cannot persist.^ This destruction of transverse angular momentum is per definition equal to a torque. Slonczewski's 
spin-transfer torque is therefore equivalent to the absorption of a spin current at an interface between a normal 
metal and a ferromagnet whose magnetization is transverse to the spin current polarization. Each electron carries an 
electric charge — e and an angular momentum of ±h/2. The loss of transverse spin angular momentum at the normal 
metal- ferromagnet interface is therefore h [I s — (I s • m) m] / (2e), where the spin-current I s is measured in the units of 
an electrical current, e.g. in Ampere. In the macrospin approximation the torque has to be shared with all magnetic 
moments or M S V of the ferromagnetic particle or film with volume V. The torque on magnetization equals the rate 
of change of the total magnetic moment of the magnet d(mM s V) stt /dt, which equals the spin current absorptiorP 
.The rate of change of the magnetization direction therefore reads: 



dm 

~dt 



2eM,V 



m x (m x I s 



(3) 



We still need to evaluate the spin current that can be generated, e.g., by the inverse spin Hall effect in the normal 
metal or optical methods. Here we concentrate on the layered normal metal-ferromagnet systems in which the current 
generated by an applied bias is polarized by a second highly coercive magnetic layer as in the schematic Fig. [T] 
Magnetoelectronic circuit theory is especially suited to handle such a problem 2 ^ For simplicity we disregard here 
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FIG. 1: Illustration of the spin-transfer torque in layered normal metaljferromagnet system, 
normal metal induces a spin-transfer torque r st t on the ferromagnet. 



A spin accumulation V^' in the 
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extrinsic dissipation of spin angular momentum due to spin-orbit coupling and disorder, which can taken into account 
when the need arises^l42l. We allow for a non-equilibrium magnetization or spin accumulation in the normal 

(s) (s) 

metal layer. V y N ' is a vector pointing in the direction of the local net magnetization, whose modulus Vjy is the 
difference between the differences in electric potentials (or electrochemical potentials divided by 2e) of both spin 
species. Including the charge accumulation vff^ (local voltage), the potential experienced by a spin-up (spin-down) 

t (c) (s) ( I (c) (s)\ 

electron along the direction of the spin accumulation in the normal metal is = -I- I Vjy = — J . 

Inside a ferromagnet, the spin accumulation must be aligned to the magnetization direction = mV F s \ Since 
V F does not directly affect the spin-transfer torque at the interface we disregard it for convenience here (see Ref. 
^] for a complete treatment), but retain the charge accumulation V F . We can now compute the torque at the 

(s) 

interface between a normal metal and a ferromagnet arising from a given spin accumulation ~V N ■ Ohm's Law for 
the spin-current projections aligned (/^) and anti-aligned (i^) to the magnetization direction then rea d 21 l 44 J (positive 
currents correspond to charge flowing from the normal metal towards the ferromagnet) 



7+ — G-f 



(4 c) -^ c) )+m.( V ^-m^ c) )], (4) 
(4 c) -l4 c) )-m.(v^-m<)]. (5) 



where and G^ are the spin-dependent interface conductances. The total charge current 1^ = , is continuous 

across the interface, iff = iff = 1^ ■ The (longitudinal) spin current defined by Eqs. Q and |i| (I-f — Jj.) m is 
polarized along the magnetization direction. The transverse part of the spin current can be written as the sum of two 

(s) 

vector components in the space spanned by the m ,V M plane as well as its normal. The total spin current on the 
normal metal side close to the interface readf^^Il 

jCa.H-) = (Jt _ h)m _ 2G W m x ( m x V W) - 2G« (m x vff) , (6) 



where G± and G^ff are two independent transverse interface conductances. blas ) j s driven by the external bias 

ViJ and should be distinguished from the pumped spin current addressed below. (R) and (I) refer to the real and 

imaginary parts of microscopic expression for these "spin mixing" interface conductances G-j-j, = G^" + iG^ff . 

The transverse components are absorbed in the ferromagnet within a very thin layer. Detailed calculations show 
that transverse spin-current absorption in the ferromagnet happens within a nanometer from the interface, where 
disorder suppresses any residual oscillations that survived the above-mentioned destructive interference in ballistic 
structures^. Spin-transfer in transition metal based multilayers is therefore an interface effect, except in ultrathin 
ferromagnetic films^. As discussed above, the divergence of the transverse spin current at the interface gives rise to 
the torque 



.(bias) _ 7^ 

stt ~ eM s V 



Gf V x (m x V# ) + G ( ff (m x V ( ff 



(7) 



Adding this torque to the Landau-Lifshitz-Gilbert equation leads to the Landau-Lifshitz-Gilbert-Slonczewski (LLGS) 
equation 

rh = —7m x H c ff + r^t™' + am x rh. (8) 

The first term in Eq. < 7 1 is the (Slonczewski) torque in the ^m, plane, which resembles the Landau-Lifshitz 

damping in Eq. (0). When the spin- accumulation is aligned with the effective magnetic field H c ff , the Slonczewski 
torque effectively enhances the damping of the ferromagnet and stabilizes the magnetization motion towards the 

(s) 

equilibrium direction. On the other hand, when is antiparallel to H ff, this torque opposes the damping. 

When exceeding a critical value it leads to precession or reversal of the magnetization. The second term in Eq. Q 

proportional to G^ff modifies the magnetic field torque and precession frequency. While the in-plane torque leads 
to dissipation of the spin accumulation, the out-of-plane torque induces a precession of the spin accumulation in the 
ferromagnetic exchange field along m. It is possible to implement the spin-transfer torque into the Landau-Lifshitz 
equation, but the conductance parameters differ from those in Eq. Q. 

Since spin currents can move magnetizations, it is natural to consider the reciprocal effect, viz. the generation of 
spin currents by magnetization motion. It was recognized in the 1970's that spin dynamics is associated with spin 
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FIG. 2: Spin-pumping in normal metal| ferromagnet systems. A dynamical magnetization "pumps" a spin current l' s ' into an 
adjacent normal metal. 

currents in normal metals. Barnes^ studied the dynamics of localized magnetic moments embedded in a conducting 
medium. He showed that the dynamic susceptibility in diffuse media is limited by the spin-diffusion length. Janossy 
and Monocf^ and Silsbee et aP^ postulated a coupling between a dynamic ferromagnetic magnetization and a spin 
accumulation in adjacent normal metals in order to explain that microwave transmission through normal metal foils is 
enhanced by a coating with a ferromagnetic la yer . The scattering theory for spin currents induced by magnetization 
dynamics was developed by Tserkovnyak et al. u on the basis of the theory of adiabatic quantum pumpinjfi^, hence 
the name "spin pumping" . Theoretical results were confirmed by the agreement of the spin-pumping induced increase 
of the Gilbert damping with experiments by Mizukami et al and Heinrich et alP^Gll. A schematic picture of spin- 
pumping in normal | ferromagnet systems is shown in Fig. [2] At not too high excitations and temperatures, the 
ferromagnetic dynamics conserves the modulus of the magnetization M s m. Conservation of angular momentum then 
implies that the spin current i^' pump ) pumped out of the ferromagnet has to be polarized perpendicularly to m, viz. 
m • j(^P u m P) _ o. Furthermore, the adiabatically pumped spin current is proportional to |rh|. Under these conditions, 
therefore, 121111 

e^.pump) = G ,(R) (m x A) + G ,{I)^ (g) 

where Gf 1 and G _[ are two transverse conductances that depend on the materials. Here the sign is defined to be 

negative when l^ ,pump ) implies loss of angular momentum for the ferromagnet. For |rh| ^ 0, the right-hand side of 
the LLGS equatio n (j8| ) must be augmented by Eq. ([9]). The leakage of angular momentum leads e.g. to an enhanced 
Gilbert damping23Ilir 

Onsager's reciprocity relations dictate that conductance parameters in thermodynamically reciprocal processes 
must be identical when properly normalized. We prove below that spin-transfer torque Q and spin pumping ^ 

indeed belong to this category and must be identical, viz. G± — G , j R and G± = G_l . Spin-transfer torque 
and spin-pumping are therefore opposite sides of the same coin, at least in the linear response regime. Since spin- 
mixing conductance parameters governing both processes are identical, an accurate measurement of one phenomenon 
is sufficient to quantify the reciprocal process. Magnetization dynamics induced by the spin-transfer torque are not 
limited to macrospin excitations and experiments are carried out at high current levels that imply heating and other 
complications. On the other hand, spin-pumping can be directly detected by the line-width broadening of FMR 
spectra of thin multilayers. In the absence of two-magnon scattering phenomena and a sufficiently strong static 
magnetic field, FMR excites only the homogeneous macrospin mode, allowing the measurement of the transverse 
conductances G'j^' and, in principle, G'jy. G'j 1 ^. Experimental results and first-principles calculations^^ agree 
quantitatively well. Rather than attempting to measure these parameters by current-induced excitation measurements, 
the values G'j_ R ^ and G 1 ^ should be inserted, concentrating on other parameters when analyzing these more complex 
magnetization phenomena. Finally we note that spin mixing conductance parameters can be derived as well from 
static magnetoresistance measurements in spin valve:P^ or by detecting the spin current directly by the inverse spin 
Hall effecPini. 
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2. Continuous Systems 



The coupling effects between (spin-polarized) electrical currents and magnetization dynamics also exist in magneti- 
zation textures of bulk metallic ferromagnets. Consider a magnetization that adiabatically varies its direction in space. 
The dominant contribution to the spin-transfer torque can be identified as a consequence of violation of angular mo- 
mentum conservation: In a metallic ferromagnet, a charge current is spin polarized along the magnetization direction 
to leading order in the texture gradients. In the bulk, i.e. separated from contacts by more than the spin-diffusion 
length, the current polarization is P = (cr-j- — a^) / (er^ + cr^), in terms of the ratio of the conductivities for majority and 
minority electrons, where we continue to measure spin currents in units of electric currents. We first disregard spin-flip 
processes that dissipate spin currents to the lattice. To zeroth order in the gradients, the spin current y a ' flowing 
is a specified (say x-) direction at position r is polarized along the local magnetization, (r) = m(r)j( s )(r). The 
gradual change of the magnetization direction corresponds to a divergence of the angular momentum of the itinerant 
electron subsystem, d x ]^ — j^d x m + md x j( s \ where the latter term is aligned with the magnetization direction 
and does not contribute to the magnetization torque. This change of spin current does not leave the electron system 
but flows into the magnetic order, thus inducing a torque on the magnetization. This process does not cause any 
dissipation and the torque is reactive, as can be seen as well from its time reversal symmetry. To first order in the 
texture gradient, or adiabatic limit, and for arbitrary current direction d 51 * 52 ! 

r( t b r ) (r) = ^^(j-V)m, (10) 

where j is the charge current density vector and the superscript "bias" indicates that the torque is induced by a voltage 



bias or electric field. From symmetry arguments another torque should exist that is normal to Eq. (10 1, but still 
perpendicular to the magnetization and proportional to the lowest order in its gradient. Such a torque is dissipative, 
since it changes sign under time reversal. For isotropic s ystem s, we can parameterize the out-of-plane torque by a 
dimensionless parameter (3 such that the total torque readiP^^, 

^tt aS) (r) - P^oP [(E • V) m + /Jmx (E ■ V) m] , (11) 

2eM s 

we have used Ohm's law, j = ctE. In the adiabatic limit, i.e. to the first order in the gradient of the magnetization 



diirij, the spin-transfer torque Eq. (11) describes how the magnetization dynamics is affected by currents in isotropic 
ferromagnets. 



Analogous to discrete systems, we may expect a process reciprocal to (11) in ferromagnetic textures similar to the 
spin pumping at interfaces. Since we are now operating in a ferromagnet, a pumped spin current is transformed into 
a charge current. To leading order a time-dependent texture is expected to pump a current proportional to the rate 
of change of the magnetization direction and the gradient of the magnetization texture. For isotropic systems, we can 
express the expected charge current as 

i| pUmp) = ^P' [m x c^m+^m] • m, (12) 

where P' is a polarization factor and /3' an out-of-plane contribution. Note that we have here been assuming a strong 
spin-flip rate so that the spin-diffusion length is much smaller than the typical length of the magnetization texture. 
Volovik considered the opposite limit of weak spin-d issip ation and kept track of currents in two independent spin 
bands^. In that regime he derived the first term in (12), proportional to P' and proved that P = P' . This results 



was re-derived by Barnes and MaekawaP^. The last term, proportional to the /3-factor was first discussed by Duine in 
Ref. (|5^|) for a mean-field model, demonstrating that (3 = f3' . More general textures and spin relaxation regimes were 
treated by Tserkovnyak and Mecklenburg^. In the following we demonstrate by the Onsager reciprocity relations 



that the coefficients appearing in the spin-transfer torques (11) are identical to those in the pumped current (12), i.e. 
P = P' and (3 = f3'. 

The proposed relations for the spin-transfer torques and pumped current in continuous systems form a local rela- 
tionship between torques, current, and electric and magnetic fields. For ballistic systems, this is not satisfied since the 
current at one spatial point depends on the electric field in the whole sample or global voltage bias and not just on the 
local electric field. The local assumption also breaks down in other circumstances. The long-range magnetic dipole 
interaction typically breaks a ferromagnet into uniform domains. The magnetization gradually changes in the region 
between the domains, the domain wall. When the domain wall width is smaller than the phase coherence length or 
the mean free path, one should replace the local approach by a global strategy for magnetization textures in which 
the dynamics is characterized by one or more dynamic (soft) collective coordinates {£ a ( T )} that are allowed to vary 
(slowly) in time 

m(rr) = m st (r;{e a (r)}), (13) 
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where m st is a static description of the texture. In order to keep the discussion simple and transparent we disregard 
thermoelectric effects, which can be important in principle^. The thermodynamic forces are —dF/d^ a , where F is 
the free energy as well as the bias voltage across the sample V. In linear response the rate of change of the dynamic 
collective coordinates and the charge current in the system are related to the thermodynamic forces —dF/d£ and V 
by a response matrix 

( \ = t in y \( - a w "i 



l)-\l H ln)\ V J' (14) 

where L^y describes the bias voltage-induced torque and Lj^ the current pumped by the moving magnetization 
texture. These expressions are general and includes e.g. effects of spin-orbit interaction. Onsager's reciprocity 
relations imply Lj^^m, H} = L^i{— m, — H} or Lj^^m^H} = L^x{— m, — H} depending on how the collective 
coordinates transform under time-reversal. The coefficient Lj^ can be easily expressed in terms of the scattering 
theory of adiabatic pumping as discussed below. This strategy was employed to demonstrate for (Ga,Mn)As that the 
spin-orbit interaction can enable a torque arising from a pure charge current bias in Ref. [42] and to compute /3 in Ref. 



3. Self- consistency: Spin-battery and enhanced Gilbert Damping 

We discussed two reciprocal effects: torque induced by charge currents (voltage or electric field) on the magnetization 
and the current induced by a time-dependent magnetization. These two effects are not independent. For instance, 
in layered systems, when the magnetization precesses, it can pump spins into adjacent normal metal. The spin- 
pumping affects magnetization dynamics depending on whether the spins return into the ferromagnet or not. When 
the adjacent normal metal is a good spin sink, this loss of angular momentum affects the magnetization dynamics 
by an enhanced Gilbert damping. In the opposite limit of little or no spin relaxation in an adjacent conductor of 
finite size, the pumped steady-state spin-current is canceled by a diffusion spin current arising from the build-up of 
spin accumulation potential in the adjacent conductor. The build-up of the spin accumulation can be interpreted as 
a spin battery^. Similarly, in magnetization textures, the dynamic magnetization pumps currents that in turn exert 
a torque on the ferromagnet. 

In the spin-battery the total spin-current in the normal metal consists of the diffusion-driven Eq. ([6]) and the pumped 
Eq. (JoJ) spin currents^. When there are no other intrinsic time-scales in the transport problem (e.g. instantaneous 
diffusion) and in the steady state, conservation of angular momentum dictates that the total spin-current in the 
normal metal must vanish, 

T (s,bias) T (s,pump) _ „ 
L N + 1 N — u > 

which from Eqs. ^ and Q results in a spin accumulation, which can be called a spin-battery bias or spin-motive 
force: 

eV^' = hm x m. (15) 

This is a manifestation of Larmor' theorem^. In diffusive systems, the diffusion of the pumped spins into the normal 
metal takes a finite amount of time. When the typical diffusion time is longer than the typical precession time, the 
AC component averages out to zerd^. In this regime, the spin-battery bias is constant and determined by 

(DC) f dt 

m x fim, (16) 



eV (s) 

N 



where r p is the precession period. Without spin-flip processes, the magnitude of the steady-state spin bias is governed 

(s) 

by FMR frequency of the magnetization precession eV y N ' = Hlofmr and is independent of the interface properties. 

Spin-flip scattering in the normal metal reduces the spin bias eV^" 1 < fiwpMR in a non-universal wa> l 15 l 56 l The loss 
of spin angular momentum implies a damping torque on the ferromagnet. Asymmetric spin- flip scattering rates in 
adjacent left and right no rmal metals can also induced a charge potential difference resulting from the spin-battery, 
which has been measured! 57 * 58 ^ The spin-battery effect has also been measured via the spin Hall effect in RefP^l. 

In the opposite regime, when spins relax much faster than their typical injection rate into the adjacent normal 
metal, ([3]), the net spin-current is well described by the spin-pumping mechanism. According to Eq. ([9|, in which 
primes may be removed because of the Onsager reciprocity, 



(pump) T^ 1 

Tstt = 2e 2 ALV 



Gf ] m x m + G i [ ) m 



(17) 



9 



We use the superscript "pump" to clarify that this torque arises from the emission of spins from the ferromagnet. 
The first term in Eq. (17) is equal to the Gilbert damping term in the LLG equation ([!]). This implies that the spin 



pumping into an adjacent conductor maximally enhances the Gilbert damping by 

(p Um p) = 7«L_ G («). (18) 
stt 2e 2 M s V ± V ' 

This damping is proportional to the interface conductance G]_ and thus the normal metal-ferromagnet surface area as 
well as inversely proportional to the volume of the ferromagnet and therefore scales as 1/dp, where dp is the thickness 
of the ferromagnetic layer. The transverse conductance per unit areas agrees well with theory^. The mi crosc opi c 



expression for G^ > and therefore q;^" 11113 '' > 0. The second term on the right hand side of Eq. (17) in (17), 

modifies the gyro-magnetic ratio and Wfmr- F° r conventional ferromagnets like Fe, Ni, and Co, G^P <C G~± by near 

cancellation of positive and negative contributions in momentum space. In these systems G^P is much smaller than 

and the effects of G^p might therefore be difficult to observe. 
A similar argument leads us to expect an enhancement of the Gilbert damping in magnetic textures. By inserting the 



pumped current Eq. ( 12 ) into the torque Eq. ( 11 ) in place of trE, we find a contribution caused by the magnetization 
dynamicfPHlU 



Ts? t rlft ^(r) = ^ P A a f([m x d i m.+j3d l m\ • m) + /3mx ([m x <9 i m+/3<9 i m] ■ rhj)l c^m, (19) 



which gives rise to additional dissipation of the order ^nrP cr/4e M S A^, where is the typical length scale for the 



variation of the magnetization texture such as the domain wall width or the radius of a vortex. Eq. ( 19 ) inserted 
into the LLG equation also renormalizes the gyromagnetic ratio by an additional factor /3. The additional dissipation 
becomes important for large gradients as in narrow domain walls and close to magnetic vortex center d 60 * 62 l 

Finally, we point out that the fluctuation-dissipation theorem dictates that equilibrium spin-current fluctuations 
associated with spin-pumping by thermal fluctuations must lead to magnetization dissipation. This connection was 
worked out in Ref. 1551 



C. Onsager Reciprocity Relations 

The Onsager reciprocity relations express fundamental symmetries in the linear response matrix relating thermo- 
dynamic forces and currents. In normal metal|ferromagnetic heterostructures, a spin accumulation in the normal 
metal in contact with a ferromagnet can exert a torque on the ferromagnet, see Eq. ([7]). The reciprocal process is 
spin pumping, a precessing ferromagnet induces a spin current in the adjacent normal metal as described by Eq. (Joj) . 
Both these effects are non-local since the spin-transfer torque on the ferromagnet arises from the spin accumulation 
potential in the normal metal and the pumped spin current in the normal metal is a result of the collective magneti- 
zation dynamics. In bulk ferromagnets, a current (or electric field) induces a spin-transfer torque on a magnetization 
texture. The reciprocal pumping effect is now an electric current (or emf) generated by the texture dynamics. In the 
next two subsections we provide technical details of the derivation of the Onsager reciprocity relations under these 
circumstance 



1. Discrete Systems 

As an example of a discrete system, we consider a normal metal-ferromagnet bilayer without any spin-orbit inter- 
action (see Ref. |3H for a more general treatment that takes spin-flip processes into account) and under isothermal 
conditions (the effects of temperature gradients are discussed in Refs. 1331641651) . The spin-transfer physics is induced 
by a pure spin accumulation in the normal metal, whose creation does not concern us here. The central ingredients 
for the Onsager's reciprocity relations are the thermodynamic variables with associated forces and currents that are 
related by a linear response matrhP^. In order to uniquely define the linear response, currents J and forces X have 
to be normalized such that F — ^XJ.. This is conventionally done by the rate of change of the free energy in the 
non-equilibrium situation in terms of currents and forces^. 

Let us consider first the electronic degrees of freedom. In the normal metal reservoir of a constant spin accumulation 
the rate of change of the free energy Fn in terms of the total spin (in units of electric charge e) reads 

F N = -s N -V ( £>. (20) 
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(s) 

This identifies V N as a thermodynamic force that induces spin currents I s = sn, which is dehned to be positive 
when leaving the normal metal. In the ferromagnet, all spins are aligned along the magnetization direction m. 

(s) 

The associated spin accumulation potential V F can only induce a contribution to the longitudinal part of the spin 
current, e.g. a contribution to the spin-current along the magnetization direction m. In our discussion of the Onsager 
reciprocity relations, we will set this potential to zero for simplicity and disregard associated change in the free energy, 
but it is straightforward to include the effects of a finite Vp^ 

Next, we address the rate of change of the free energy related to the magnetic degrees of freedom in the ferromagnet, 

F(m) = -M s VH cS • m/T, 

where F(m) is the magnetic free energy. The total magnetic moment M s Vm is a thermodynamic quantity and the 
effective magnetic field H c ff = —dF/d(M s Vm) is the thermodynamic force that drives the magnetization dynamics 
m. 

In linear response, the spin current I s = s and magnetization dynamics M s Vm. are related to the thermodynamic 
forces as 



M s Vm. 



fj^ram) Jpns) 



(21) 



where L {mm \ U ms \ U sm \ and are 3x3 tensors in, e.g., a Cartesian basis for the spin and magnetic moment 
vectors. Onsager discovered that microscopic time-reversal symmetry leads to relations between the off-diagonal 
components of these linear-response coefficients. Both magnetization in the ferromagnet and the spin-accumulation 
in the normal metal are anti-symmetric under time-reversal leading to the reciprocity relations 



L{r>(m) = 4"»>(-m). 



(22) 



Some care should be taken when identifying the Onsager symmetries in spin accumulation-induced magnetization 
dynamics. Specifically, the LLCS equation ^ cannot simply be combined with the linear response relation (21 1 and 
Eq. (22). Only the Landau-Lifshitz-Slonczewski (LL) Eq. ^ directly relates m to H e g as required by Eq. (21 1. In 
terms of the 3x3 matrix O e.g. 



Oij(m) 



k 



(23) 



where ey/- = \ (j — i) {k — i) (k — j) is the Levi-Civita tensor, m x H G g = OH ff , and the LLGS |8| equation can be 



written as 



(l-aO)m = 0(- 7 H eff ) + T stt . 



(24) 



By Eq. (21), the pumped current in the absence of a spin accumulation (Vj^ = 0) is iff = 
9|), I^y = X^^rh, where the 3x3 tensor x( sm ) has components 



Xl sm) (m) = -- 
e 



'j 



<->"_|_ / J e inj m n T <~r± / J £ink m n£klj m k 



nkl 



(25) 



From the LLG equation (24) for a vanishing spin accumulation (vff = 0) and thus no bias-induced spin-transfer 



torque (r. 



bias) 



0), the pumped spin current can be expressed as iff = X^ sm ^O 



identifies the linear response coefficient Z,( sm ) in terms of X^ sm ^ as 



1 - aO 



(-7 H eff), which 



= - 7 x (sm ^d \l-aO 



(26) 



Using the Onsager relation (22) and noticing that Oy(m) = Oj*(— m) and X^ m '(m) — JQ^ mJ (— m) 



= _ 7 I -aO OX^ 



-l _ 



(27) 
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The rate of change of the magnetization by the spin accumulation therefore becomes 



_J_£(m.) v W 



N ' 



7 



1 - «o 



OX 



(sm)y(s) 



N 



(28) 



.( drift ) 



Furthermore, the LLGS equation (24) in the absence of an external magnetic field reads 1 — aO 
Inserting the phenomenological expression for the spin-transfer torque ([7]), we identify the linear response coefficient 



.(drift) 
stt 



7 AyMyH 

M,V N • 



7 



M,Ve 



Gfmx 



m x V 



This agrees with the phenomenological expression ^ when 



Q'iR) 



r {R). 



G 



'(^) 



m x V 



(29) 



(30) 



Spin-pumping as expressed by Eq. ([9| is thus reciprocal to the spin-transfer torque as described by Eq. Q. In 
Sec jUI A 1| these relations are derived by first principles from quantum mechanical scattering theory, resulting in. 
q'(R) _ _ ( e 2 /h) J2 nm ^nm — r nm ( r nm)* f° r a narrow constriction, where r\ m (r^ m ) is the reflection coefficient 
for spin-up (spin-down) electrons from waveguide m to waveguide mode n . For layered systems with a constant cross 
section the microscopic expressions of the transverse (mixin g) conductances should be renormalized by taking into 
account the contributions from the Sharvin resistance d 23 * 66 !, which increases the cond uctan ce by roughly a factor of 
two and is important for a quantitatively comparison between theory and experiments!^^ 



2. Continuous Systems 



The Onsager reciprocity relations also relate the magnetization torques and currents in the magnetization texture 
of bulk magnets. Following Refs. (|31I32[) . the rate of change of the free energy related to the electronic freedom in the 
ferromagnet is Fp = — J drqV, where q is the charge density and eV = fi is the chemical potential. Inserting charge 
conservation, q + V • j = and by partial integration 



F F = 



drj-E 



(31) 



which identifies charge as a thermodynamic variable, while the electric field E = S7V is a thermodynamic force which 
drives the current density j. For the magnetic degrees of freedom, the rate of change of the free energy (or entropy) is 

(32) 



F m = -M s l drm(r)-H efI (r). 

Just like for discrete systems, H c ff (r), is the thermodynamic force and Mgm is the thermodynamic variable to which 
it couples. In a local approximation the (linear) response depends only on the force at the same location: 



M s m 
j 



jjjnra) fftnE) 
l(Em) l(EE) 



A/ s H c ff 
E 



(33) 



where U mm \ U m i\ L^ m \ and £C») 

are the local response functions. Onsager's reciprocity relations dictate again 

that 



:H m) (m) = L^(-m). 



(34) 



Starting from the expression for current pumping (12), we can determine the linear response coefficient L^ Em ^ from 



= -J-^-crP' [ejkim k dimi+f3 r diirij] , 
ij Ze 



(35) 
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FIG. 3: Schematic of how transport between a normal metal and a ferromagnet is computed by scattering theory. The 
scattering region, which may contain the normal metal-ferromagnet interface and diffusive parts of the normal metal as well 
as ferromagnet, is attached to real or fictious leads that are in contact with a left and right reservoir. In the reservoirs, the 
distributions of charges and spins are assumed to be known via the charge potential and spin accumulation bias. 



where the operator O is introduced in the same way as for discrete systems (23) to transform the LLG equation into 
the LL form p3|. According to Eq. ([34]) 



O" 1 



1 - aO 



= -IttvP' [eikimkdjmi-P'djmi 
ij Ze 



(36) 



The change in the magnetization induced by an electric field is then M s rh 



(bias) 



£,( m J')E so that the spin-transfer 



torque due to a drift current r 



(bias) 



1 - aO 



(bias) 



can be written as 



.(bi: 



' 2eM. 



aP' 

f-trnn 



r£f S) = 7§Jf [(E • V) m + /3'm x E ■ Vm] 



(37) 
(38) 



This result agrees with the phenomenological expression for the pumped current (12) when P = P' and f3 = /?'. 
Therefore, the pumped current and the spin-transfer torque in continuous systems are reciprocal processes. The 
pumped current can be formulated as the response to a spin-motive forced. 

In small systems and thin wires, the current-voltage relation is not well represented by a local approximation. 
A global approach based on collective coordinates as outlined around Eq. ( 13 ) is then a good choice to keep the 
computational effort in check. Of course, the Onsager reciprocity relations between the pumped current and the 
effective current-induced torques on the magnetization hold then as welP^l. 



III. MICROSCOPIC DERIVATIONS 



A. Spin-transfer Torque 



1. Discrete Systems - Magneto-electronic Circuit Theory 



Physical properties across a scattering region can be expressed in terms of the region's scattering matrix, which 
requires a separation of the system into reservoirs, leads, and a scattering region, see Fig. Q. In the lead with index 
a, the field operator for spin s-electrons is^ 



f de 




1 


- v (ns) 



-1/2 



no 



>t/h 



(e) + e~ ikx b^ s \e) 



(39) 



in terms of the annihilation operators a^ s ' (fri™^) for particles incident on (outgoing from) the scattering region 
in transverse wave guide modes with orbital quantum number n and spin quantum number s (s =f or s =1). 
Furthermore, the transverse wave function is ip^^^g), the transverse coordinate g, the longitudinal coordinate along 
the waveguide is x and vi™^ is the longitudinal velocity for waveguide mode ns. The positive definite momentum k is 



related to the energy e by hk — (2me) 1 ^ 2 . 
by the scattering matrix 



The annihilation operators for incident and outgoing electrons are related 



/3ms' 



(eK ; (e) 



(40) 
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In the basis of the leads (a = N (normal metal) or a = F (ferromagnet)), the scattering matrix is 

where r (t) is a matrix of the reflection (transmission) coefficients between the wave guide modes for an electron 
incident from the left. Similarly, r and t characterize processes where the electron is incident from the right. 



In terms of the field operators defined by Eq. (39) and the scattering matrix Eq. (40), at low frequencies, the spin 
current that flows in the normal metal a = N in the direction towards the scattering region is 

/DC pOO 
dei / ^EEE^'^-^^A^^'^^^a^^eOaCK)^), (41) 
-OO J — OO a l i 

P7 nml era' 



where 



j^(nm.nl)((T,a / ) ^ 



»( eil£2 ) = [<W* (nm) * (sCT) ^ 7 <5 ( *^ (sV) - S% s - ma) *{e 1 )S^'^'\e 2 ) 

ss' 

and cr^ ss - 1 is a vector of the 2x2 Pauli matrices that depends on the spin indices s and s' of the waveguide mode. 
The charge current can be found in a similar way. We are interested in the expectation value of the spin-current (41 ) 
when the system is driven out-of-cquilibrium. In equilibrium, the expectation values are 

(aL ns)t (e)4 mS 'V)) oq = S(e-e')S aP d (ss '^^f FB (e), (42) 

where /fd(c) is the Fermi-Dirac distribution of electrons with energy e. A non-equilibrium spin-accumulation in the 



normal metal reservoir is not captured by the local equilibrium ansatz in Eq. (42), however. A spin accumulation in 
the normal metal reservoir can still be postulated when spin-flip dissipation is slow compared to all other relevant time 
scales. We assume the normal metal and ferromagnet have an isotropic distribution of spins in the orbital space, and 
for clarity consider no charge bias. The expectation for the number of charges and spins in the waveguide describing 
normal metal leads attached to the normal reservoirs are 

(a^e^V)) = 5(e - e') [^">^')/ PD (e) + 5^f^' s \e)] . (43) 
The spin-accumulation is related to the 2x2 out-of-equilibrium distribution matrix s \e) by 

poo 

«t("').V«=/ de$'\e)/e. (44) 



For the spin-transfer physics, a bias voltage in the ferromagnet does not contribute since it only gives rise to a charge 
current and a longitudinal spin current. As in the previous section, we therefore set this voltage to zero for simplicity, 
so that in the ferromagnetic lead attached to the ferromagnetic reservoir 

(4" s)t ( £ )4 ms V)) = S(e - e> (ms) <5 (s ' s) /FD( £ ). (45) 

Furthermore, the expectation values of the cross-correlations remain zero also out-of-equilibrium, 
a^ ts ^(e)a^ ls (e')/ = 0- The spin current in lead a is then 



/OO 
-OO I I I 



(46) 



nmlss' era' 

Without spin-flip scattering, the reflection coefficient can be expressed as 

r" = (r^ T + rT/) 5^/2 + m • a sa (fg$ - r™/) /2 (47) 
which can be represented in spin space as 

~NN — r NN 1 + N N m " 
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since the scattering matrix can be decomposed into components aligned and anti-aligned with the magnetization 
direction. These matrices only depend on the orbital quantum nu mb ers (n and to). Using the representation of the 
out-of-equilibrium spin density in terms of the spin accumulation (44 1 21 , 



I# = (G t + G;)m(m.V 



^G^m x 



m x V 



2G 



2>( 



m x V 



(49) 



in agreement with (6]) when there is no bias voltage in the ferromagnet (Vf = 0) which we have assumed for clarity 
here. We identify the microscopic expressions for the conductances- 2 — associated with spins aligned and anti-aligned 
with the magnetization direction 



Gt 



Gt 



E 

nm 

E 



nm,t 
NN 



nm,t 
NN 



and the transverse (complex valued) spin-mixing conductance 



G± 



El 



c nm^ nm,J.* 

Onm r NN r NN 



(50) 
(51) 

(52) 



These results are valid when the transmission coefficients are small such that currents do not affect the reservoirs. 
Otherwise, the transverse cond uctance parameters should be renormalized by taking into account the Sharvin resis- 
tances, as described above^EH. In the limit we considered here, the expression for the spin-current depends only on 
the reflection coefficients for transport from the normal metal towards the ferromagnet and not on the transmission 
coefficients for propagation from the normal metal into the ferromagnet. This follows from our assumption that the 
ferromagnet is longer than the transverse coherence length as w ell as our disregard of the spin accumulation in the 
ferromagnet. Both assumptions can be easily relaxed if necessar y ^ 15 * 19 *. 



2. Continuous Systems 

Spin torques in continuous spin textures can be studied by either quantum kinetic theory,^ imaginary-time^ 
and functional KeldystF^ diagrammatic approaches, or the scattering-matrix formalism.^ The latter is particularly 
powerful when dealing with nontrivial band structures with strong spin-orbit interactions, while the others give com- 
plementary insight, but are mostly limited to simple model studies. When the magnetic texture is sufficiently smooth 
on the relevant length scales (the transverse spin coherence length and, in special cases, the spin-orbit precession 
length) the spin torque can be expanded in terms of the local magnetization and current density as well as their 
spatial-temporal derivatives. An example is the phenomenological Eq. ( |11[ ) for the electric-field driven magnetization 
dynamics of an isotropic ferromagnet. While the physical meaning of the coefficients is clear, the microscopic origin 
and magnitude of the dimensionless parameter /3 has still to be clarified. 

The solution of the LLG equation ([lj appended by these spin torques depends sensitively on the relationship 
between the dimensionless Gilbert damping constant a and the dissipative spin-torque parameter f3: the special case 
P/a = l effectively manifests Galilean invariance^ while the limits ^/a>l and j3/a <C 1 are regimes of qualitatively 
distinct macroscopic behavior. The ratio (3/ a determines the onset of the ferromagnetic current-driven instability^ 
as well as the Walker threshold^ for the current-driven domain- wall motiorP^, and both diverge as (3 /a — > 1. The 
sub-threshold current-driven domain- wall velocity is proportional to /3/ap3 while f3/a — 1 in a special point, at which 
the effect of a uniform current density j on the magnetization dynamics is eliminated in the frame of reference that 
moves with velocity v cx j, which is of the order of the electron drift velocity.^ Although the exact ratio f3/a is a 
system-dependent quantity, some qualitative aspect s not too sensitive to the microscopic origin of these parameters 
have been discussed in relation to metallic systems! 68 * 69 * 71 * 74 -! However, these approaches fail for strongly spin-orbit 
coupled systems such as dilute magnetic semiconductors^!. 

Let us outline the microscopic origin of (3 for a simple toy model for a ferromagnet. In Ref. 1681 we developed a 
self-consistent mean-field approach, in which itinerant electrons are described by a single-particle Hamiltonian 

U = [Uo + U(r, t)} 1 + ^<x • (H + H xc ) (r, t) + U a , (53) 

where the unit matrix 1 and a vector of Pauli matrices <x = (a x , a y , a z ) form a basis for the Hamiltonian in spin 
space. "Hp is the crystal Hamiltonian including kinetic and potential energy. U is the scalar potential consisting of 
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disorder and applied electric-field contributions. The total magnetic field consists of the applied, H, and exchange, 
H xc , fields that, like U , are parametrically time dependent. Finally, the last term in the Hamiltonian, % a , accounts 
for spin-dephasing processes, e.g, due to quenched magnetic disorder or spin-orbit scattering associated with impurity 
potentials. This last term is responsible for low-frequency dissipative processes affecting dimensionless parameters a 
and /3 in the collective equation of motion. 

In the time-dependent spin-density-functional theorjPsHZU of itinerant ferromagnetism, the exchange field H xc is a 
functional of the time-dependent spin-density matrix 

p afl (r,i) = <*J(r)* a (r)) t , (54) 

where ^>'s are electronic field operators, which should be computed self-consistently as solutions of the Schrodingcr 
equation for %. The spin density of conducting electrons is given by 

s(r) = ^Tr[«rp(r)] . (55) 

We focus on low-energy magnetic fluctuations that are long ranged and transverse and restrict our attention to a 
single parabolic band. Consideration of more realistic band structures is also in principle possible from this starting 
poinlPS. We adopt the adiabatic local-density approximation (ALDA, essentially the Stoner model) for the exchange 
field: 

7 HH xc [p}{r,t) w A xc m(r,t), (56) 



with direction m = — s/s locked to the time-dependent spin density (55). 

In another simple model of ferromagnetism, the so-called s-d model, conducting s electrons interact with the 
exchange field of the d electrons that are assumed to be localized to the crystal lattice sites. The d-orbital electron 
spins account for most of the magnetic moment. Because d-electron shells have large net spins and strong ferromagnetic 
correlations, they are usually treated classically. In a mean-field s-d description, therefore, conducting s orbitals are 



described by the same Hamiltonian ( 53 ) with an exchange field ( 56 ) . The differences between the Stoner and s-d 
models for the magnetization dynamics are subtle and rather minor. In the ALDA/Stoner model, the exchange 
potential is (on the scale of the magnetization dynamics) instantaneously aligned with the total magnetization. In 
contrast, the direction of the unit vector m in the s-d model corresponds to the d magnetization, which is allowed 
to be slightly misaligned with the s magnetization, transferring angular momentum between the s and d magnetic 
moments. Since most of the magnetization is carried by the latter, the external field H couples mainly to the d 



spins, while the s spins respond to and follow the time-dependent exchange field (56). As A xc is usually much larger 
than the external (including demagnetization and anisotropy) fields that drive collective magnetization dynamics, the 
total magnetic moment will always be very close to m. A more important difference of the philosophy behind the 
two models is the presumed shielding of the d orbitals from external disorder. The reduced coupling with dissipative 
degrees of freedom would imply that their dynamics are more coherent. Consequently, the magnetization damping has 
to originate from the disorder experienced by the itinerant s electrons. As in the case of the itinerant ferromagnets, 
the susceptibility has to be calculated self-consistently with the magnetization dynamics parametrized by m. For 
more details on this model, we refer to Refs. I7§1 and 1551 With the above differences in mind, the following discussion 
is applicable to both models. The Stoner model is more appropriate for transition-metal ferromagnets because of 
the strong hybridization between d and s,p electrons. For dilute magnetic semiconductors with by deep magnetic 
impurity states the s-d model appears to be a better choice. 

The single-particle itinerant electron response to electric and magnetic fields in Hamiltonian ( |53[ ) is all that is 
needed to compute the magnetization dynamics microscopically. Stoner and s-d models have to be distinguished only 
at the final stages of the calculation, when we self-consistently relate m(r, t) to the electron spin response. The final 



result for the simplest parabolic-band Stoner model with isotropic spin-flip disorder comes down to the torque (11) 
with a rs /3. The latter is proportional to the spin-dephasing rate r" 1 of the itinerant electrons: 

(57) 

The derivation assumes uj,^ 1 <C A xc /fi, which is typically the case in real materials sufficiently below the Curie 



temperature. The s-d model yields the same result for /3, Eq. (57), but the Gilbert damping constant 

a rs 77/3 (58) 
is reduced by the ratio 77 of the itinerant to the total angular momentum when the d-electron spin dynamics is not 



damped. [Note that Eq. (58) is also valid for the Stoner model since then 77 = 1.] 
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These simple model considerations shed light on the microscopic origins of dissipation in metallic ferromagnet as 
reflected in the a and j3 parameters. In Sec. IV we present a more systematic, first-principle approach based on the 
scattering-matrix approach, which accesses the material dependence of both a and f3 with realistic electronic band 
structures. 



B. Spin Pumping 



Discrete Systems 



When the scattering matrix is time-dependent , the energy of outgoing and incoming states does not have to be 
conserved and the scattering relation (40) needs to be appropriately generalized^. We will demonstrate here how 



this is done in the limit of slow magnetization dynamics, i.e., adiabatic pumping. When the time dependence of the 

a,(nrr 
'a/3 



scattering matrix S^ l J n ' > [Xi(t)] is parameterized by a set of real- valued parameters Xi(t), the pumped spin current in 



excess of its static bias-driven value ( 49 ) is given bjB^l 

i 

where the "spin emissivity" vector by the scatterer into lead a is^ 



(59) 



dX r 



dS. 



ms.ncr )* 



( 

a/3 



13 mn ss'cr 



dX, 



_£.(ss') Q(ms',ncr) 



a/3 



(60) 



Here, <j( ss ' is again the vector of Pauli matrices. In the case of a magnetic monodomain insertion and in the absence 
of spin-orbit interactions, the spin-dependent scattering matrix between the normal-metal leads can be written in 
terms of the respective spin-up and spin-down scattering matrices!^ 



(61) 



Here, m(i) is the unit vector along the magnetization direction and f (i) are spin orientations defined along (opposite) 
to m. 



Spin pumping due to magnetization dynamics m(t) is then found by substituting Eq. (61) into Eqs. (60) and ([59 
After straightforward algebraP^ 



Kit) 



- 1 dt x dt J 



(62) 



As before, we assume here a sufficiently thick ferromagnet, on the scale of the transverse spin-coherence length. Note 
that the spin pumping is expressed in terms of the same complex-valued mixing conductance G± = G 



iG { [ ] as 



the dc current ( 49 1 , in agreement with the Onsager reciprocity principle as found on phenomenological grounds in 
Sec. HTCl 



Charge pumping is governed by expressions similar to Eqs. (59) and (60), subject to the following substitution: 



<t — y S (Kronecker delta). A finite charge pumping by a monodomain magnetization dynamics into normal-metal leads, 
however, requires a ferromagnetic analyzer or finite spin-orbit interactions and appropriately reduced symmetries, as 
discussed in Refs. SHHH 



An immediate consequence of the pumped spin current ( 62 ) is an enhanced Gilbert damping of the magnetization 
dynamics.^ Indeed, when the reservoirs are good spin sinks and spin backflow can be disregarded, the spin torque 
associated with the spin current (62 ) into the a-th lead, as dictated by the conservation of the spin angular momentum, 
Eq. contributes (cf. Eq. |18j)): 



2e 2 M S V 



(63) 



to the Gilbert damping of the ferromagnet in Eq. ([T]). Here, g* ~ 2 is the g factor of the ferromagnet, M S V its total 

magnetic moment, and hb is Bohr magneton. For simplicity, we neglected G^[\ which is usually not important for 
inter-metallic interfaces. If we disregard energy relaxation processes inside the ferromagnet, which would drain the 
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associated energy dissipation out of the electronic system, the enhanced energy dissipation associated with the Gilbert 
damping is associated with heat flows into the reservoirs. Phenomenologically, the dissipation power follows from the 
magnetic free energy F and the LLG Eq. ([!]) as 



P = -d m F m ■ m = A/ s VH cff • rh = 



aM s V 

7 



m 



or, more generally, for anisotropic damping (with, for simplicity, an isotropic gyromagnetic ratio), by 

F = m am. 



7 



(64) 



(65) 



Heat flows can be also calculated microscopically by the scattering-matrix transport formalism. At low tempera- 
tures, the heat pumping rate into the a-th lead is given bj^ 



jE _ |^(ms,ns') 
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(66) 



where the carets denote scattering matrices with suppressed transverse-channel indices. When t he time dependence 

Mm 
7 ap 



is entirely due to the magnetization dynamics, g^ s ^ ns ) — 



dmS (rns,ns) . m uti l izing 

again Eq. fl61b, we find for the 



heat current into the a-th leadP^ 
in terms of the dissipation tensor^ 



I a = m • G a ■ m , 



-in — ' y <';//, orrij 



(67) 



(68) 



In the limit of vanishing spin-flip in the ferromagnet, meaning that all dissipation takes place in the reservoirs, we 
find 
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Equating this 1^ with P above, we obtain a microscopic expression for the Gilbert damping tensor a: 



a 



2e 2 M,V ' 



(69) 



(70) 



which agrees with Eq. (63). Indeed, in the absence of spin-orbit coupling the damping is necessarily isotropic. While 
Eq. (63) reproduces the additional Gilbert damping due to the interfacial spin pumping, Eq. (69) i s mor e general, 
and can be used to compute bulk magnetization damping, as long as it is of a purely electronic origi n 88 ^ 89 !. 



2. Continuous Systems 

As has already been noted, s pin pumping in continuous systems is the Onsager counterpart of the spin-transfer 



torque discussed in Sec. Ill A 2 2D While a direct diagrammatic calculation for this pumping is possible^, with re- 
sults equivalent to those of the quantum-kinetic description of the spin-transfer torque outlined above, we believe 
that the scattering-matrix formalism is the most powerful microscopic approach^. The latter is particularly suit- 
able for implementing parameter-free computational schemes that allow a realistic description of material-dependent 
properties. 

An important example is pumping by a moving domain wall in a quasi-one-dimensional ferromagnetic wire. When 
the domain wall is driven by a weak magnetic field, its shape remains to a good approximation unaffected, and only 
its position r w (t) along the wire is needed to parameterize its slow dynamics. The electric current pumped by the 
sliding domain wall into the a-th lead can then be viewed as pumping by the r w parameter, which leads to 81 
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The total heat flow into both leads induced by this dynamics is according to Eq. ( 66 1 



dr w dr w 




(72) 



Evaluating the scattering-matrix expressions on the right-hand side of the above equations leads to microscopic 
magnetotransport response coefficients that describe the interaction of the domain wall with electric currents, including 
spin transfer and pumping effects. 

These results leads to microscopic expressions for the phenomenological response^! of the domain-wall velocity f w 
and charge current I c to a voltage V and magnetic field applied along the wire H: 



2 AM S H 
V 



(73) 



subject to appropriate conventions for the signs of voltage and magnetic field and assuming a head-to-head or tail- 
to-tail wall such that the magnetization outside of the wall region is collinear with the wire axis. 2AM S H is the 
thermodynamic force normalized to the entropy production by the magnetic system, where A is the cross-sectional 
area of the wire. We may therefore expect the Onsager's symmetry relation L cw — L wc . When a magnetic field 
moves the domain wall in the absence of a voltage I c = (L cw / L ww )r w , which, according to Eq. (71 ) leads to the ratio 
L cw /L ww in terms of the scattering matrices. The total energy dissipation for the same process is I E = r^ /L ww , 



which, according to Eq. (72), establishes a scattering-matrix expression for L ww alone. By supplementing these 



equations with the standard Landauer-Biittiker formula for the conductance 



G — ^Tr (^S\ 2 Si2^j , 

valid in the absence of domain- wall dynamics, we find L cc in the same spirit since G = L c 
the phenomenological response coefficients in Eq. ( 73 1 reacP^l 



dr w dr„ 
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(74) 

L 2 wc / L ww Summarizing, 

(75) 

(76) 
(77) 



When the wall is sufficiently smooth, we can model spin torques and pumping by the continuum theory based on 
the gradient expansion in the magnetic texture, Eqs. (Ill and (12). Sol ving f or the magnetic- field and current-driven 
dynamics of such domain walls is then possible using the Walker ansatz^Eni. Introducing the domain- wall width A„, : 



2AMM, 



and j3 = 



hPG L v 



(78) 



When the wall is sharp the adiabatic approximation underlying the leading-order gradient expansion breaks down. 
These relations can still be used as definitions of the effective domain- wall a and j3. As such, these could be distinct 
from the bulk values that are associated with smooth textures. This is relevant for dilute magnetic semiconductors, 
for which the adiabatic approximation easily breaks dowrP^. In transition-metal ferromagnets, on the other hand, 
the adiabatic approximation is generally perceived to be a good starting point, and we may expect the dissipative 
parameters in Eq. (78) to be comparable to their bulk values discussed in Sec. Ill A 2 



IV. FIRST-PRINCIPLES CALCULATIONS 



We have shown that the essence of spin pumping and spin transfer can be captured by a small number of phe- 
nomenological parameters. In this section we address the material dependence of these phenomena in terms of the 
(reflection) mixing conductance G±, the dimensionless Gilbert damping parameter a, and the out-of-plane torque 
parameter /3. 
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For discrete systems the (reflection) mixing conductance G± was studied theoretically by Xia et alW^, Zwierzycki et 
aiP^ and Carva et alW^. G± describes the spin current flowing in response to an externally applied spin accumulation 
eV s that is a vector with length equal to half of the spin-splitting of the c hemical p otentials e|V s | = e(V^ — Vj,)/2. It 
also describes the spin torque exerted on the moment of the magnetic layeJ 9 ^^ 45 * 91 * ^. Consider a spin accumulation 
in a normal metal N, which is in contact with a ferromagnet on the right magnetized along the z axis. The spin 
current incident on the interface is proportional to the number of incident channels in the left lead, lf n = 2G^ h V s , 
while the reflected spin current is given by 



I N - 2 

A out — z 



G { [ ] Gf-G^ 1 1 | V.. ("") 

V o o 




where G a are the conventional Landauer-Biittiker conductances. The real and imaginary parts of G^ h — G± = 
(e 2 /h) J2mn r mn r mn are related to the components of the reflected transverse spin current and can be calculated by 
considering a single N|F interface^. When the ferromagnet is a layer with finite thickness d sandwiched between 
normal metals, the reflection mixing conductance depends on d and it is necessary to consider also the transmission 
mixing conductance (e 2 /h) Yl mn Knn^'nm- ^ n Ref- H5J both reflection and transmission mixing conductances were 
calculated for Cu|Co|Cu and Au|Fe|Au sandwiches as a function of magnetic layer thickness d. The real and imaginary 
parts of the transmission mixing conductance and the imaginary part of the reflection mixing conductance were shown 
to decay rapidly with increasing d implying that the absorption of the transverse component of the spin current occurs 
within a few monolayers of the N|F interface for ideal lattice matched interfaces. When a minimal amount of interface 
disorder was introduced the absorption increased. The limit G± — > GjJ/ 1 corresponds to the situation where all of the 
incoming transverse polarized spin current is absorbed in the magnetic layer. The torque is then proportional to the 
Sharvin conductance of the normal metal. This turns out to be the situation for all but the thinnest (few monolayers) 
and cleanest Co and Fe magnetic layers considered by Zwierzycki et aiP^ However, when there is nesting between 
Fermi surface sheets for majority and minority spins so that both spins have the same velocities over a large region of 
reciprocal space, then the transverse component of the spin current does not damp so rapidly and Gj_ can continue 
to oscillate for large values of d. This has been found to occur for ferromagnetic Ni in the (001) direction.^ 

Eq. [TJJ implies that the spin pumping renormalizes both the Gilbert damping parameter a and the gyromagnetic 
ratio 7 of a ferromagnetic film embedded in a conducting non-magnetic medium. However, in view of the results 
discussed in the previous paragraph, we conclude that the main effect of the spin pumping is to enhance the Gilbert 
damping. The correction is directly proportional to the real part of the reflection mixing conductance and is essentially 
an interface property. Oscillatory effects are averaged out for realistic band structures, especially in the presence of 

{ R) 

disorder. G]_ determines the damping enhancement of a single ferromagnetic film embedded in a perfect spin-sink 
medium and is usually very close to G^ h for intermetallic interface d 91 * 93 !. 



A. Alpha 

We begin with a discussion of the small-angle damping measured as a function of temperature using ferromagnetic 
resonance (FMR). There is general agreement that spin-orbit coupling and disorder are essential ingredients in any 
description of how spin excitations relax to the ground state. In the absence of intrinsic disorder, one might expect the 
damping to increase monotonically with temperature in clean magnetic materials and indeed, this is what is observed 
for Fe. Heinrich et a/P^l developed an explicit model for this high-temperature behaviour in which itinerant s electrons 
scatter from localized d moments and transfer spin angular momentum to the lattice via spin-orbit interaction. This s— 
d model results in a damping that is inversely proportional to the electronic relaxation time, a ~ 1/r, i.e., is resistivity- 
like. However, at low temperatures, both Co and Ni exhibit a sharp rise in damping as the temperature decreases. 
The so-called breathing Fermi surface model was proposecPsHSHl to describe this low-temperature conductivity-like 
damping, a ~ t. In this model the electronic population lags behind the instantaneous equilibrium distribution due 
to the precessing magnetization and requires dissipation of energy and angular momentum to bring the system back 
to equilibrium. 

Of the numerous microscopic models that h ave been proposed^ to explain the damping behaviour of metals, only 
the so-called "torque correlation model" (TCM^P^is qualitatively successful in explaining the non-monotonic damping 
observed for hep Co that results from conductivity-like and resistivity-like behaviours at low and high temperatures, 
respectively. The central result of the TCM is the expression 
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for the damping. The commutator [a-,"H S o] describes a torque between the spin and orbital moments that arises 



as the spins precess. The corresponding matrix elements in (80) describe transitions between states in bands n 
and m induced by this torque whereby the crystal momentum k is conserved. Disorder enters in the form of a 
phcnomenological relaxation time r via the spectral overlap 

W»,m(k) = -- [ A n (e,k)A m (e,k)^fde (81) 

7T J CLE 

where the electron spectral function A n (e,\s) is a Lorentzian centred on the band n, whose width is determined by 
the scattering rate. For intraband transitions with m = n, integration over energy yields a spectral overlap which is 
proportional to the relaxation time, like the conductivity. For interband transitions with m ^ n, the energy integration 
leads to a spectral overlap that is roughly inversely proportio nal to th e relaxation time, like the resistivity. 

To interpret results obtained with the TCM, Gilmore et al.E^SHUHl used an effective field approach expressing the 
effective field about which the magnetization precesses in terms of the total energy 

«H- . - (82) 

and then approximated the total energy by a sum of single particle eigenvalues E ~ k £„k/nk, so that the effective 
field naturally splits into two parts 



H cff = 1 
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(83) 



the first of which corresponds to the breathing Fermi surface model, intraband transitions and conductivity-like 
behaviour while the second term could be related to interband transitions and resistivity-like behaviour. Evaluation 
of this model for Fe, Co and Ni using first-principles calculations to determine £ n k including spin-orbit coupling 
yields result s for the damping a in good qualitative and reasonable quantitative agreement with the experimental 
observations 

In spite of this real progress, the TCM has disadvantages. As currently formulated, the model can only be applied 
to periodic lattices. Extending it to handle inhomogeneous systems such as ferromagnetic substitutional alloys like 
Permalloy (NisoFe2o), magnetic multilayers or heterojunctions, disordered materials or materials with surfaces is far 
from trivial. The TCM incorporates disorder in terms of a relaxation time parameter r and so suffers from the 
same disadvantages as all transport theories similarly formulated, namely, that it is difficult to relate microscopically 
measured disorder unambiguously to a given value of t. Indeed, since r in general depends on incoming and scattered 
band index n, wave vector k, as well as spin index, assuming a single value for it is a gross simplification. A useful 
theoretical framework should allow us to study not only crystalline materials such as the ferromagnetic metals Fe, Co 
and Ni and substitutional disordered alloys such as permalloy (Py), but also amorphous materials and configurations 
such as magnetic heterojunctions, multilayers, thin films etc. which become more important and are more commonly 
encountered as devices are made smaller. 

The scattering theoretical framework discussed in section III B satis fies these requirements and has recently been 
implemented by extending a first-principles scattering formalis m 1 105 ! 106 ! based upon the local spin density approxima- 
tion (LSDA) of density functional theory (DFT) to include non-collinearity, spin-orbit coupling (SOC) and chemical 
or thermal disorder on equal footings'®' Relativistic effects are inc luded by using the Pauli Hamiltonian. To calculate 
the scattering matrix, a "wave-function matching" (WFM ) scheme^SHIoZI implemented with a minimal bas is of tight- 
binding linearized muffin-tin orbitals (TB-LMTOsJP2SEnil Atomic-sphere-approximation (A SA) potential d 108 * 109 ! are 
calculated self-consistently using a surface Green's function (SGF) method also implemented^! with TB-LMTOs. 



1. NiFe alloys. 



The flexibility of the scattering theoretical formulation of transport can be demonstrated with an application to 
NiFe binary alloys! 8 -^ Charge and spin densities for binary alloy A and B sites are calculated using the coherent 
potential approximation (CPA) n -i< generalized to layer structure*^. For the transmission matrix calculation, the 
resulting spherical potentials are dist ributed at random in large lateral supercells (SC) subject to maintenance of the 
appropriate concentration of the allo} * 105 * 1 - ^. Solving the transport problem using lateral supercells makes it possible 
to go beyond effective medium approximations such as the CPA. As long as one is only interested in the properties 
of bulk alloys, the leads can be chosen for convenience and Cu leads with a single scattering state for each value of 
crystal momentum, kii are very convenient. The alloy lattice constants are determined using Vegard's law and the 
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FIG. 4: Calculated resistivity as a function of the concentration x for fee Nii-^Fe^ bina ry alloy s with (solid line) and without 
(dashed-dotted line) SOC. Low temperature experimental results are shown as symbolP^EZI. The composition NigoFe2o is 
indicated by a vertical dashed line. Inset: resistance of Cu|NigoFe2o|Cu as a function of the thickness of the alloy layer. Dots 
indicate the calculated values averaged over five configurations while the solid line is a linear fit. 



lattice constants of the leads are made to match. Though NiFe is fee only for the concentration range < x < 0.6, 
the fee structure is used for all values of x. 

To illustrate the methodology, we begin by calculating the electrical resistivity of NigoFe2o- In the Landa uer-Biit tiker 
formalism, the conductance can be expressed in terms of the transmission matrix t as G — (e 2 /h)Tr {ttf J1U2U12J. The 
resistance of the complete system consisting of ideal leads sandwiching a layer of ferromagnetic alloy of thickness L 
is R(L) = l/G(L) = l/Crsh + 2-Rif + Rb(L) where Gsh = (2e 2 //i) N is the Sharvin conductance of each lead with N 
conductance channels per spin, i?;t is t he inte rface resistance of a single N|F interface, and Rh{L) is the bulk resistance 
of a ferromagnetic layer of thickness l) 66 l 106 l. When the ferromagnetic slab is sufficiently thick, Ohmic behaviour is 
recovered whereby Rh(L) ~ pL as shown in the inset to Fig. [4] and the bulk resistivity p can be extracted from the 
slope of R(L). For currents parallel and perpendicular to the magnetization direction, the resistivities are different 
and have to be calculated separately. The average resistivity is given by p = (py + 2pj_)/3, and the anisotropic 
magnetoresistance ratio (AMR) by (p\\ — p±)/p. 

For NigoFe2o we find values of p = 3.5 ± 0.15 /iOhm-cm and AMR = 19 ± 1% , compared to experimental low- 
temperature values in the range 4.2 — 4.8 /iOhm-cm for p a nd 18% for AMRP^. The resistivity calculated as a 
function of x is co mpared to low temperature literature values^^^ in Fig. [4] The overall agreement with previous 
calculations is g 00c JllS1119| j n S pite of the smallness of the SOC, the resistivity of Py is underestimated by more than 
a factor of four when it is omitted, underlining its importance for understanding transport properties. 

Assuming that the Gilbert damping is isotropic for cu bic substitutional alloys and allowing for the enhancement 
of the damping due to the F|N interface d 14 * 17 * 18 * 45 * 120 * 121 !, the total damping in the system with a ferromagnetic slab 
of thickness L can be written G(L) = Git + Gb(L) where we express the bulk damping in terms of the dimensionless 
Gilbert damping parameter Gb(L) = ar)M s (L) — a"//i s AL, where p s is the magnetization density and A is the cross 
section. The results of calculations for Ni 8 oFe2o arc shown in the inset to Fig. [5) The intercept at L = 0, Gtf, 
allows us to extract the damping enhancement^ but here we focus on the bulk properties and leave consideration 
of the material dependence of the interface enhancement for later study. The value of a determined from the slope 
of G(L)/ (~fp s A) is 0.0046 ± 0. 0001 that is at the lower end of the range of values 0.004 — 0.013 measured at room 
temperature for P y 17 ! 18 * 121 ™ . 

Fig. [5] shows the Gilbert damping parameter as a function of x for Nii^Fe^ binary alloys in the fee structure. 
From a large value for clean Ni, it decreases rapidly to a minimum at x ~ 0.65 and then grows again as the limit 
of clean fee Fe is approached. Part of the decrease in a with increasing x can be explained by the increase in the 
magnetic moment per atom as we progress from Ni to Fe. The large values of a cal culated in the dilute alloy limits 
can be understood in terms of conductivi ty- like enhancement at low temperature s] 132 * ^ 3 -! that has been explained 
in terms of intraband scatteringj^QSJHlQSlISl] xhe trend exhibited by the theoretical a(x) is seen to be reflected by 
experimental results obtained at room temperature. In spite of a large spread in measured values, these seem to be 
systematica lly larg er than the calculated values. Part of this discrepancy can be attributed to an increase in a with 
temperatur e 122 * 134 !. 

Calculating a for the end members, Ni and Fe, of the substitutional alloy Nii-^Fea, presents a practical problem. In 
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FIG. 5: Calculated zero temperature (solid line) and experimental room t emperature (symbols) values of the Gilbert damping 
parameter as a function of the concentration x for fee Nii^Fe^ binary alloys^- 18 * 120 * ^^. Inset: total damping of Cu|Ni8oFe2o|Cu 
as a function of the thickness of the alloy layer. Dots indicate the calculated values averaged over five configurations while the 
solid line is a linear fit. 
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FIG. 6: Calculated Gilbert damping and resistivity for fee Ni as function of the relative RMS displacement with respect to the 
corresponding lattice constant, ao = 3.524 A. 



these limits there is no scattering whereas in experiment there will always be some residual disorder at low temperatures 
and at finite temperatures, electrons will scatter from the thermally displaced ions. We introduce a simple "frozen 
thermal disorder" scheme to study Ni and Fe and simulate the effect of temperature via electron-phonon coupling by 
using a random Gaussian distribution of ionic displacements u^, corresponding to a harmonic approximation. This 
is characterized by the root-mean-square (RMS) displacement A = yj (|uj| 2 ) where the index i runs over all atoms. 
Typical values will be of the order of a few hundredths of an angstrom. We will not attempt to relate A to a real 
lattice temperature here. 

We calculate the total resistance R(L) and Gilbert damping G(L) for thermally disordered scattering regions of 
variable length L and extract the resistivity p and damping a from the slopes as before. The results for Ni are 
shown as a function of the RMS displacement in Fig. (6) The resistivity is seen to increase monotonically with A 
underlining the corr elation between A and a real temperature. For large val ues of A, a saturates for Ni in agreement 
with experiment!^] an( j calculations based on the torque-correlation modet ^ * 103 * 104 ! where no concrete scattering 
mechanism is attached to the relaxation time r. The absolute value of the saturated a is about 70% of the observed 
value. For small values of A, the Gilbert damping increases rapidly as A decreases. This sharp rise corresponds to the 
experimentally observed conductivity-like behaviour at low temperatures and confirms that the scattering formalism 
can reproduce this feature. 
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FIG. 7: (a) Sketch of the configuration of a Neel DW in Py sandwiched by two Cu leads. The arrows denotes local magnetization 
directions. The curve shows the mutual angle between the local magnetization and the transport direction [z axis), (b) 
Magnetization profile of the rotated Neel wall, (c) Magnetization profile of the Bloch wall. 



B. Beta 



To evaluate expressions ( 78 ) for the out-of-plane spin-torque parameter j3 given in Section IIIB requires modelling 



domain walls (DW) in the scattering region sandwiched between ideal Cu leads. A head-to-head Neel DW is introduced 
inside the permalloy region by rotating the local magnetization to follow the Walker profile, Jn(z) = [f(z),0,g(z)] 
with f(z) = cosh - [(2 — r w )/X w ] and g(z) = — tanh[(z — r w )/X w ] as shown schematically in Fig.r7[a). r w is the DW 
center and A^, is a parameter characterizing its width. In addition to the Neel wall, we also study a rotated Neel wall 
with magnetization profile m(z) = [g(z), 0, f(z)] sketched in Fig. ^b) and a Bloch wall with m(z) = [g(z), f(z), 0] 
sketched in Fig. [7jc). 




FIG. 8: Calculated effective Gilbert damping constant a for Py DWs as a function of A„. The dashed lines show the calculated 
a for bulk Py with the magnetization parallel to the transport direction®'. 
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FIG. 9: Calculated out-of-plane spin torque parameter Pf3 for permalloy DWs as a function of A^, 



The effective Gilbert damping constant a of permalloy in the presence of all three DWs calculated using ( 78 ) is 
shown in Fig. [8] For different types of DWs, a is identical wit hin the numerical accuracy indicating that the Gilbert 
damping is isotropic due to the strong impurity scattering^®!. In the adiabatic limit, a saturates to the same value 
(the dashed lines in Fig. [8]) calculated for bulk permalloy using ( [68] ) . It implies that the DWs in permalloy have 
little effect on the magnetization relaxation and the strong impurity scattering is the dominant mechanism to release 
energy and magnetization. This is in contrast to DWs in (Ga,Mn)As where Gilbert damping is mostly contributed by 
the reflection of the carriers from the DWP^ At A^ < 5 nm, the non-adiabatic reflection of conduction electrons due 
to the rapidly-varying magnetization direction becomes significant and results in a sharp rise in a for narrow DWs. 

The out-of-plane torque is formulated as (3(h-fP/2eM s )m x (j • V)m in the Landau-Lifshitz-Gilbert (LLG) equation 
under a finite current density j. In principle, the current polarization P is required to determine f3. Since the spin- 
dependent conductivities of permalloy depend on the angle between the current and the magnetization, P is not 
well-defined for magnetic textures. Instead, we calculate the quantity P/3, as shown in Fig. [9] for a Bloch DW. For 
X w < 5 nm, P/3 decreases quite strongly with increasing A^ corresponding to an expected non-adiabatic contribution 
to the out-o f-pla ne torque. This arises from the spin-flip scattering induced by the rapidly-varying magnetization in 
narrow DWs^ 5 -' and does not depend on the specific t ype of DW. For A^, > 5 nm, which one expects to be in the 
adiabatic limit, P/3 decreases slowly to a constant valu d 27 * 32 * 53 * 69 * 78 *^M143|. It is unclear what length scale is varying 
so slowly. Unfortunately, the spread of values for different configurations is quite large for the last data point and our 
best estimate of P/3 for a Bloch DW in permalloy is ~ 0.08. Taking the theoretical value of P ~ 0.7 for permallojP^ 
our best estimate of (3 is a value of ~ 0.01. 



V. THEORY VERSUS EXPERIMENTS 



Spin-torque induced magnetization dynamics in multilayers and its reciprocal effect, the spin pumping, are ex- 
perimentally well established an d qua ntitatively understood within the framework described in this paper, and need 
not to be discussed further here |l£ 122l R ecent FMR experiments also confirm the spin-pu mping contribution to the 
enhanced magnetization dissipation^ 4 ! Spin-pumping occurs in magnetic insulators as wel l 7 * 145 ! 

The parameters that control the current-induced dynamics of continuous textures are much less well known. Most 
experiments are carried out on permalloy (Py) . It is a magnetically very soft material with large domain wall widths 
of the order of 100 nm. Although the adiabatic approximations appears to be a safe assumption in Py, many systems 
involve vortex domain walls with large gradients in the wall center, and, therefore, possibly sizable nonadiabatic 
corrections. Effective description for such vortex dynamics has been constructed in Ref. I62|, where it was shown, 
in particular, that self-consistent quadratic corrections to damping (which stem from self-pumped currents inducing 
backaction on the magnetic order) i s generally non-negligible in transit ion- metal ferromagnets. 

Early experimental studies^ 46 * 147 ! for the torque-supplemented [Eq. (11)] LLG equation describing current-driven 
domain- wall motion in magnetic wires reported values of the j3/a ratios in Py close to unity, in agreement with simple 
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Stoner- model calculations. However, much larger values /3/a ~ 8 was extracted from the current-induced oscillatory 
motion of domain walls The inequality/3 7^ a was also inferred from a cha ract eristic tran sverse to vortex wall 
structure transformation, although no exact value of the ratio was established. 14 ^ 1 ' In RefP^I, vanadium doping of 
Py was shown to enhance j3 up to nearly 10a, with little effect on a itself. Even larger ratios, /3/a ~ 20, were 
found for magnet ic vortex motion by an analysis of their displacement as a function of an applied dc current in disc 
structures.™^ 

Eltschka et o/P^ reported on a measurement of the dissipative spin- torque parameter f3 entering Eq. (11), as 
manifested by a thermally-activated motion of transverse and vortex domain walls in Py. They found the ratio 
0v/@t ~ 7 for the vortex vs transverse wall, attributing the larger f3 to high magnetization gradients in the vortex wall 
core. Their ratio fit 1 01 ~ 1.3 turns out to be close to unity, where a is the bulk Gilbert damping. The importance of 
large spin-texture gradients on the domain-wall and vortex dynamics was theoretically discussed in Refs. 1601621 

The material dependence of the current-induced torques is not yet well investigated. A recent study on CoNi and 
FePt wir es w ith perpendicular magnetization found j3 ~ a, in spite of the relatively n arro w domain walls in these 
materials P^ Current-induced domain-wall dynamics in dilute magnetic semiconductors^^ generally exhibit similar 
phenomenology, but a detailed discussion, especially of the domain wall creep regime that can be accessed in these 
systems, is beyond the scope of this review. 

Finally, the first term in the spin-pumping expression ( fl2| ) has been measured by Yang et oiP^ for a domain wall 
moved by an applied magnetic field above the Walker breakdown field. These experiments confirmed the existence 
of pumping effects in magnetic textures, which are Onsager reciprocal of spin torques and thus expected on general 
grounds. Similar experiments carried out below the Walker breakdown would also give direct access to the f3 parameter. 



VI. CONCLUSIONS 



A spin polarized current can excite magnetization dynamics in ferromagnets via spin-transfer torques. The reciprocal 
phenomena is spin-pumping where a dynamic magnetization pumps spins into adjacent conductors. We have discussed 
how spin-transfer torques and spin-pumping are directly related by Onsager reciprocity relations. 

In layered normal metal-ferromagnet systems, spin-transfer torques can be expressed in terms of two conductance 
parameters governing the flow of spins transverse to the magnetization direction and the spin-accumulation in the nor- 
mal metal. In metallic systems, the field-like torque is typically much smaller than the effective energy gain/damping 
torque, but in tunnel systems they might become comparable. Spin-pumping is controlled by the same transverse 
conductance parameters as spin-transfer torques, the magnetization direction and its rate of change. It can lead to 
an enhanced magnetization dissipation in ultra-thin ferromagnets or a build-up of spins, a spin-battery, in normal 
metals where the spin-flip relaxation rate is low. 

Spin-transfer torque and spin-pumping phenomena in magnetization textures are similar to their counterparts in 
layered normal metal-ferromagnet systems. A current becomes spin polarized in a ferromagnet and this spin-polarized 
current in a magnetization texture gives rise to a reactive torque and a dissipative torque in the lowest gradient 
expansion. The reciprocal pumping phenomena can be viewed as an electromotive force, the dynamic magnetization 
texture pumps a spin-current that in turn is converted to a charge current or voltage by the giant magnetoresistance 
effect. Naturally, the parameters governing the spin-transfer torques and the pumping phenomena are also the same 
in continuously textured ferromagnets. 

When the spin-orbit interaction becomes sufficiently strong, additional effects arise in the coupling between the 
magnetization and itinerant electrical currents. A charge potential can then by itself induce a torque on the ferro- 
magnet and the reciprocal phenomena is that a precessing ferromagnet can induce a charge current in the adjacent 
media. The latter can be an alternative way to carry out FMR measurements on small ferromagnets by measuring 
the induced voltage across a normal metal-ferromagnet-normal metal device. 

These phenomena are well-know and we have reviewed them in a unified physical picture and discussed the con- 
nection between these and some experimental results. 
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